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Quantum trajectory tests of radical-pair quantum dynamics in CIDNP measurements 

of photosynthetic reaction centers 
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Chemically induced dynamic nuclear polarization is a ubiquitous phenomenon in photosynthetic 
reaction centers. The relevant nuclear spin observables are a direct manifestation of the radical-pair 
mechanism. We here use quantum trajectories to describe the time evolution of radical-pairs, and 
compare their prediction of nuclear spin observables to the one derived from the radical-pair master 
equation. While our approach provides a consistent description, we unravel a major inconsistency 
within the conventional theory, thus challenging the theoretical interpretation of numerous CIDNP 
experiments sensitive to radical-pair reaction kinetics. 


I. INTRODUCTION 

Quantum coherence in the light-harvesting process of 
photosynthesis (ll-[l5l| is a central theme in the growing 
field of quantum biology [T^ . Equally important for 
understanding photosynthesis and potentially designing 
biomimetic devices harvesting solar energy is the charge 
separation following the trapping of the exciton’s energy 
in the photosynthetic reaction center. Charge separa¬ 
tion proceeds via a cascade of electron transfers between 
radical-ion pairs, as for example bacteriochlorophyll and 
bacteriopheophytin. In parallel with charge, spin also 
plays a ubiquitous role [TJ , with the effect of chemically 
induced dynamic nuclear polarization (CIDNP) mani¬ 
fested in a wide range of photosynthetic reaction centers 

[Hill- 

Nuclear spin effects in CIDNP are regulated by the 
radical-pair mechanism, which describes a class of spin- 
dependent chemical reactions studied W smn chemistry 
[Hiq . We have recently addressed the funda¬ 

mental quantum dynamics of the radical-pair mechanism, 
using concepts from quantum information science to de¬ 
scribe the intertwined effects of coherent spin motion and 
spin-dependent charge recombination of radical-pairs, ar¬ 
riving at what we understand is a fundamental master 
equation describing the time evolution of p, the radical- 
pair’s spin density matrix. The traditional (also called 
Haberkorn) theory used until now is a limiting case 
of our theory valid in the regime of strong spin relaxation. 

Interestingly, in CIDNP measurements the lifetimes of 
singlet and triplet radical-pairs (RPs) are small enough to 
allow observation of quantum effects without them being 
masked by spin relaxation. Hence CIDNP appears to be 
an ideal setting to test our understanding of the quantum 
dynamics of the radical-pair mechanism. To this end, 
we demonstrated a notable difference between our 
theory and Haberkorn’s approach in predicting a CIDNP 
effect at earth’s field. However, measurements at earth’s 
field are still challenging. Moreover, in we assumed 
equal singlet and triplet recombination rates, ks = fcx, 
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whereas in real reaction centers we encounter asymmetric 
recombination rates, kr ^ fcs [^ . 

Instead of attempting a comparison of the two differ¬ 
ent approaches based on an absolute quantitative predic¬ 
tion, we here use CINDP as a testbed to study the in¬ 
ternal consistency of the two theories, Haberkorn’s and 
ours, and we do so using the realistic (i.e. asymmetric) 
recombination rates and high magnetic fields pertinent 
to CIDNP experiments. It is a basic fact of the the¬ 
ory of open quantum systems that the time evolution 
of an ensemble of systems described by a master equa¬ 
tion must exactly reproduce the average calculated from 
single-system realizations, called quantum trajectories. 
We here use a central CIDNP observable, the nuclear 
spin polarization of the radical-pair’s reaction products, 
and show that Haberkorn’s theory produces vastly differ¬ 
ent predictions when using Haberkorn’s master equation 
compared to averaging Haberkorn’s quantum trajecto¬ 
ries, while our theory produces largely consistent predic¬ 
tions either way. 


II. QUANTUM MEASUREMENT APPROACH 

TO RADICAL-PAIR QUANTUM DYNAMICS 

A. Master equation approach 

In recent years we have addressed the fundamental 
quantum dynamics of radical-pair reactions, schemati¬ 
cally depicted in FiglTJ This biochemical system is both 
a leaky and an open quantum system, losing popula¬ 
tion due to the spin-dependent recombination reactions, 
while simultaneously suffering decoherence. The latter 
aspect of the dynamics was addressed in [s^ , and from a 
slightly different perspective in [H | , using quantum mea¬ 
surement theory. The former was phenomenologically 
considered in IdOll . while a first-principles approach was 
developed in [44| . where we formally defined a singlet- 
triplet coherence measure, called Pcoh- We also utilized 
the quantum-communications concept of quantum retro- 
diction to derive the spin-dependent reaction terms of the 
master equation. We will briefly recapitulate the results 
of [i^ for completeness of this work. 

The full master equation we arrived at in (43| , the con- 
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Open and Leaky Quantum System 



DA 


FIG. 1: A charge transfer from a photoexcited donor-acceptor 
dyad DA produces a radical-ion pair in the singlet state 
®D*’''A*~. Intramolecule magnetic interactions (Hamiltonian 
%) induce a coherent interconversion between singlet and 
triplet radical-pairs, ^ ''’D*^A*“, terminated by 

the recombination event, that results into a singlet, DA, or a 
triplet neutral product, "'^DA. Radical-pairs are leaky systems 
because real transitions to the singlet and triplet vibrational 
reservoir states lead to population loss at the recombination 
rate fcs and fex, respectively. They are also open quantum sys¬ 
tems, because virtual transitions to the same reservoir states 
and back to the radical-pair lead to singlet-triplet dephasing 
at the rate (fcs -I- fcT)/2. 


sistency of which we will explore in this work, reads 


dp 

dt 


fcs + fcx 


(pQs + Qsp — SQspQs) 


( 1 ) 

( 2 ) 


— (I — Pcoh)(fcsQspQs + fcTQTpQx) (3) 

*= + ‘'’''^(Qs.Qs + QtpQt+ 

-QspQt H-Qt/'QsI (4) 

Pcoh Vrnh ' 


Pcoh ■ 


dt 


The density matrix p describes the spin state of the 
radical-ion-pair, consisting of the two unpaired electrons 
and the nuclear spins residing in the two radicals. The 
density matrix and all other operators are represented 
by d X d matrices, where d is the dimension of the total 
spin space of the radical-pair. In this work we will con¬ 
sider a radical-pair having just one nuclear spin, hence 
d = 8, and the 8 basis kets are |S) 0 l-fl), |S) 0 IJJ-), 
|T+) 0 1^), |T+) ® 1^), |To) ® 1^), |To) ® 1^), |T_) 0 1^), 
|T_) 0 |JJ.). The two-electron state is on the left of the 
tensor product and the nuclear spin state on the right. 
The two-electron spin subspace is spanned by the sin¬ 
glet |S) = (Iti) - \\A))/\/2 and the triplets |T+) = |tt), 

|To) = (in) + in))/V2, |T_) = in). 

The operators Qs and Qt project the radical-pair spin 
state onto the electron singlet and triplet subspace, re¬ 
spectively. They are orthogonal, complete and idempo- 
tent, i.e. QsQt = QtQs = 0, Qs + Qt = 1, where 
1 is the d X d unit matrix, and Qg = Qs, Qx = Qt- 
Writing and for the spin operator of the donor’s 
and acceptor’s unpaired electron, respectively, it is Qs = 


il- Sd • SA- 

The rates fcs and are the singlet and triplet recom¬ 
bination rates. If at t = 0 we prepare an RP ensemble 
in the singlet (triplet) electron spin state, and assume 
that there is no singlet-triplet (S-T) mixing, then the RP 
population would decay exponentially at a rate fcs (fcx)- 
These two rates are properties of the particular RP un¬ 
der consideration. In principle they can be calculated 
from electron transfer theory, but in practice they are 
determined from experiment. There are RPs for which 
fcs = fcx, and there are RPs for which fcs fcx- 

The term o in the previous master equation is the 
ordinary unitary evolution driven by the intramolecule 
magnetic interactions contained in the Hamiltonian % 
(Zeeman, hyperfine etc). The particular Hamiltonian 
used in this work is relevant to CINDP measurements and 
will be outlined in Section IV. Since singlet and triplet 
states are not eigenstates of %, the term o generates 
S-T coherence. 

This is dissipated by the Lindblad term which 
we derived in |39l . 1^ . and which describes a continu¬ 
ous quantum measurement of Qs performed by the vi¬ 
brational reservoirs of FigHJ the results of which are un¬ 
observed. This measurement effects projections to the 
electron singlet or triplet subspace suffered by individual 
RPs at random times. The interruptions of the coher¬ 
ent S-T mixing at the single-molecule level lead to S-T 
decoherence at the ensemble level. 

The final two terms ([3]) and (jH) are the so-called reac¬ 
tion terms, reducing the RP population, given by Tr{p}, 
in a spin-dependent way. The normalization of the ini¬ 
tial density matrix is Trjpo} = 1- If we choose a small 
enough dt so that fcsdt, fcxdf "C I, the fraction of the RP 
population that will recombine into singlet and triplet 
neutral products within the interval dt around time t 
is drs = fcsdtTr{pQs} and drx = fcxdlTrjpQx}, re¬ 
spectively. Based on how coherent is the RP ensemble 
at time t, quantified by Proii, which is a function of p 
straightforward to calculate [ij, we use the theory of 
quantum retrodiction to probabilistically estimate the 
pre-recombination state of the observed reaction prod¬ 
ucts and thus arrive at m and ®. As expected, it is 
dTr{p} = —drs — drx. In the single-molecule picture, 
drg (drx) is the probability that a single radical-pair will 
recombine in the singlet (triplet) channel during the time 
interval dt. 

The values of pcoh range from 0, describing a maxi¬ 
mally incoherent mixture of singlet and triplet RPs, to 
0 < Pcoh < 1 describing partially coherent, to Pcoh = I 
describing maximally coherent RPs. When pcoh = 0, the 
whole term o vanishes. 

Given (i) the two rates fcs and fcx, (h) the specific 
Hamiltonian TL, and (hi) the initial state po, Eq. 0-0 
can be used to propagate p in time and thus calculate the 
time-dependence of any physical observable of interest. 
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B. Quantum trajectory approach 

The new physics of RP quantum dynamics that we in¬ 
troduced in is that during its lifetime, a radical- 

pair undergoes random projections to the singlet or 
triplet electron spin subspace. These projections take 
place at different and random times for each radical-pair. 
They run simultaneously with and independently of the 
second kind of random event, the RP charge recombina¬ 
tion, which terminates the reaction. To generate quan¬ 
tum trajectories for an RP being in the state |'0) at time 
t, we thus have to consider in total 5 possible events that 
can take place in the following time interval dt: 


TABLE I: Quantum trajectory time evolution according to 
Kominis’ approach. 


Name 

Event 

Probability of event 

KI 

projection to the singlet 

state 

v/ <'0IQsl’/’> 

dps = 

K2 

projection to the triplet 

state 

\/ <'0IQtIV'> 

dpT - 

K3 

singlet recombination 

drs = ksdt{ip\Qs\‘ip) 

K4 

triplet recombination 

dn = kTdt{tp\QT\'tp) 

K5 

hamiltonian evolution 

1 — dps — dpT — drs — drT 


It is well known from the theoretical treatment of open 
quantum systems that the master equation describ¬ 
ing the time evolution of the system’s density matrix 
should exactly reproduce the average of many single¬ 
system quantum trajectories. Thus, the average of many 
trajectories formed by K1-K5 should exactly reproduce 
the results of the master equation O-O- We will check 
whether this is the case in the context of CIDNP observ¬ 
ables in Section IV. 


III. HABERKORN APPROACH TO 
RADICAL-PAIR QUANTUM DYNAMICS 

A. Master equation approach 

The traditional, or Haberkorn master equation reads 

^ = -i[H,p] - y(QsP + pQs) - -^(QtP + pQt) (5) 

Interestingly, Haberkorn’s master equation follows from 
our master equation 0-0 by forcing pcoh to be zero at 
all times. In any case, as we have done for our approach 
in the previous section, Haberkorn’s approach must be 
able to provide the equivalent picture of single-molecule 
quantum trajectories. However, the concept of quantum 
trajectories has not been utilized in spin chemistry so far. 
Furthermore, although a recent experiment (4^ provided 
evidence for the physical reality of the S-T decoherence 


process we introduced, a general consensus on what ex¬ 
actly is the quantum state evolution of surviving RPs is 
still missing. From our perspective, we do not see how 
Haberkorn’s approach, being phenomenological, can in¬ 
corporate quantum trajectories without going into the 
derivations of [H, which lead to events KI and K2, 
but we leave it as an open question to be addressed by 
the proponents of the conventional theory. Nevertheless, 
we will here outline some rather strong guidelines as to 
how the consistency check of Haberkorn’s approach can 
in principle unfold. 


B. Quantum trajectory approach 

From Haberkorn’s theory point of view, it is clear that 
one cannot agree with possibilities KI and K2, since these 
lead to our approach. Hence one has to suggest what is 
the specific state evolution of surviving RPs, i.e. what is, 
if any, the state change of radical-pairs until the instant 
of their recombination into a neutral product. 

We will now show that there is limited freedom in do¬ 
ing so. This is because in order to secure consistency 
in the dynamically simple case ks = ^t, one has to ac¬ 
cept what has been until recently the intuitive answer 
to the previous question, namely that nothing (besides 
Hamiltonian evolution) happens to surviving RPs. As 
mentioned in Section HA, the rates ks and fcx are pa¬ 
rameters entering into the master equation, which obvi¬ 
ously must be valid for any choice of those parameters. 
The case ks = fcx is rather simple dynamically, since in 
this case RP population decays exponentially at a rate 
fc = fcs = fcx, without the decay affecting the state of the 
surviving RPs. This means the following. Consider for 
example a 50/50 mixture of singlet and triplet RPs, hav¬ 
ing no magnetic interactions (fH — 0). If ks = k^, the 
same number of singlet and triplet RPs will recombine 
in the interval dt, hence at time t + dt the mixture will 
still be 50/50, albeit having a smaller total population. 
On the other hand, if ks k^, the spin character of this 
mixture would change, becoming more (less) singlet if 
ks < fcx {ks > kr). 

To summarize, (i) assuming that having a different fun¬ 
damental theory for different radical pairs (i.e. different 
combinations of ks and fcx) is not an acceptable option, 
and (ii) being unable to propose what happens in the gen¬ 
eral case to non-recombining RPs from Haberkorn’s point 
of view, except in the special case ks = k^, where consis¬ 
tency forces one to accept that nothing else happens be¬ 
sides unitary evolution (to be proved in the following), we 
take this to be the general answer, and hence Haberkorn’s 
quantum trajectories are formed by the three events pre¬ 
sented in Table H. 

There are two comments to be made. First, the recom¬ 
bination probabilities drs and drx are the same in both 
Haberkorn’s approach and ours, since both theories agree 
in how the singlet and triplet reaction yields, Fj = f drj, 
with j = S, T, are calculated. Second, we can now easily 
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TABLE II: Quantum trajectory time evolution according to 
Haberkorn’ approach. 


Name 

Event 

Probability of event 

HI 

singlet recombination 

drs = ksdt{'p\Q,s\ip) 

H2 

triplet recombination 

drx = A:Tdt(i/) Qxh/’} 

H3 

hamiltonian evolution 

1 — drs — drT 


prove our previous statement about HI-H3 ensuring con¬ 
sistency in the special case ks = A:t- Indeed, taking into 
account the completeness relation Qs + Qt = 1, it follows 
that when ks = kT = k, Haberkorn’s master equation ([5|) 
becomes dp/dt = —i['H,p] — kp. Defining p = e~^*R, 
it follows that dR/dt = —It is thus evident 
that apart from an exponential decay of Tr{p}, the only 
radical-pair state change is due to the Hamiltonian evo¬ 
lution. In terms of the quantum trajectories H1-H3, we 
can retrieve the master equation dp/dt = —i/H, p] — kp as 
follows. Since Qs + Qt = 1, it is drs + drx = fcdtTrjp}. 
The single-RP density matrix is p/Tr{p}, hence averag¬ 
ing H1-H3 leads to 


p + dp = drsip- 
+ dvT{p — 


Tr{p} 


Tr{p} 

+ {1 - drs - drT){p - ii'HjPldt), (6) 


from which it easily follows that indeed dp/dt = 
—i\/H^p\ — kp. To reiterate, the quantum trajectories 
formed by H1-H3 exactly reproduce the master equation 
([S]) in the special case ks = kr- In the following we will 
show that this consistency check fails in the general and 
more interesting (in terms of realistic applications) case 
ks ^ kT- 


IV. TESTING THE CONSISTENCY OF 
KOMINIS’ AND HABERKORN’S APPROACHES 
USING CIDNP OBSERVABLES 

We will now demonstrate that while our approach is 
largely (but still not perfectly) consistent, Haberkorn’s 
approach is highly inconsistent. We will use a simple one- 
nuclear-spin radical-ion-pair with parameters relevant to 
CIDNP experiments. We stress that how many nuclear 
spins we consider, or which particular Hamiltonian we 
pick to exhibit the aforementioned inconsistency is of no 
concern, since as well known, it takes many supporting 
cases to establish a theory, but just one counterexample 
to invalidate it. Nevertheless, we choose a Hamiltonian 
of the same form considered in CINDP works like [l^ , 


Larmor frequency, and A and B isotromc and anisotropic 
hyperfine coupling constants. As in |26l |. we take a mag¬ 
netic field of 5 T along the z-axis. For the rest of 
the parameters we use Ag = 4 x 10“'*, A = Aw and 
B = A/2. Finally, we use the asymmetric recombination 
rates ks = (20 ns)“* and fcx = (1 ns)“* pertinent to 
photosynthetic reaction centers |17l |. 

We now calculate the nuclear spin deposited to the re¬ 
action products, diz, as a function of time [s^. 

(a) Density matrix propagation If p is the RP den¬ 
sity matrix at time t, during dt there will be drs = 
ksdtTrlpQs} singlet and dr^ = fcTdtTrlpQx} triplet 
neutral products, the properly normalized density ma¬ 
trix of which is ps = QspQs/Tr{pQs} and pt = 
QTpQT/TrlpQx}, respectively. Hence the total (singlet 
-I- triplet) ground-state nuclear spin accumulated during 
dt is 

diz = drsTr{Rps} + drxTrjRpt} (8) 

= dt(^fcsTr{RQspQs} + kxTr{RQx/oQx}) (9) 

(b) Quantum trajectories If \il}) is the state of the 
radical-pair at the random instant of recombination, 
then the properly normalized state of the singlet and 
triplet reaction product is |^s) = Qs|V’)/\/(V’lQslV’) and 
l^t) = Qt|V’)/\/('0IQt|'0)i respectively. Hence the nu¬ 
clear spin deposited to the reaction product is {ips\Iz\'>/s) 
for a trajectory terminating with a singlet recombination, 
and {ipt\Iz\'4’t) for a trajectory terminating with a triplet 
recombination. 

We evolve the quantum trajectories, and we numer¬ 
ically solve the master equation for two different initial 
conditions, (II) \ipo) = |S)ig)|JJ.), and (12) \ipo) = |S)®|D-). 
We then average the results of (II) and (12) in order to 
simulate the realistic scenario of starting with an unpo¬ 
larized nuclear spin. Finally, we integrate the averaged 
traces in order to find f dIz, which is directly accessible 
in CIDNP experiments. This is normalized by the Boltz¬ 
mann equilibrium nuclear spin, which for the considered 
parameters is /thermal ~ 10“®. 

In Figl2]we show the main result of this work. The level 
of inconsistency of Haberkorn’s approach is evident just 
by a visual inspection of Figsl2Kd-f). Even more impres¬ 
sive is the result for f dR, reproduced for convenience 
in Table HI. The results of Haberkorn’s master equation 
and Haberkorn’s quantum trajectories are not only dif¬ 
ferent in magnitude by 740%, but are also of different 
sign. In our approach the two results are of the same 
sign and different in magnitude by 50%. In Appendix A 
we describe in detail how the quantum trajectories are 
produced, while in Appendix B we elaborate on the ac¬ 
curacy of our calculations. 


Aw Aw 

R = ~;^SAz - ^SDz + aillz + As Aziz + BsAzIx, (7) 


V. DISCUSSION 


where Aw is the difference in the Larmor frequencies of The consistency of a theory is not proof of correct- 
donor and acceptor electrons due to Ag, w/ the nuclear ness, but on the contrary, the inconsistency of a theory 
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FIG. 2: Monte Carlo (grey) and master equation (black) calculation of the nuclear spin diz deposited to the neutral reaction 
products as a function of time, normalized by the thermal nuclear spin at the field of 5 T and 233 °K. (a)-(c) Kominis’ and 
(d)-(f) Haberkorn’s approach. We use two different initial conditions, |^o} = |S} (8) |11} (a,d) and lipo) = |S) ® Ifl) (b,e). In 
(c) and (f) we plot the average [(a)-|-(b)]/2 and [(d)-|-(e)]/2, respectively, in order to simulate an initially unpolarized nuclear 
spin. In the insets of (c) and (f) we show the integral of the respective traces, reflecting the total nuclear spin of the reaction 
products at the end of the reaction. 


TABLE III: Internal consistency of Haberkorn’s and Kominis’ 
approaches to radical-pair quantum dynamics tested by the 
reaction products’ integrated nuclear spin polarization rele¬ 
vant to CIDNP measurements. 


J diz j -^thermal 

Kominis 

Haberkorn 

Master Equation 

-3603 

-4250 

Monte Carlo 

-2410 

505 


is proof of its inadequacy. Hence while this work pro¬ 
vides a supporting argument that our approach is in the 
right direction, it unravels a major inconsistency within 
the traditional approach to radical-pair dynamics. 

Until now, all the interesting information about the 
physical properties of photosynthetic reaction centers 
were extracted from CIDNP signals based on the tra¬ 
ditional understanding of the radical-pair mechanism. In 
other words, several mechanisms so far understood to 


produce enhanced nuclear spin polarizations are mostly 
based on the combined action of very specific Hamilto¬ 
nian interactions and radical-pair reactions kinetics. It is 
clear that the extracted physical information contained 
in the former will be skewed by the inconsistent descrip¬ 
tion of the latter. 

In the conclusions of it was stated that ’’Until an 
experimental instance is found that requires an alterna¬ 
tive description of the recombination kinetics, we rec¬ 
ommend continued use of the conventional approach”. 
We believe that a failed consistency check of the con¬ 
ventional theory is quite stronger than an experimental 
instance challenging the theory. Even more so in light of 
the previous comment, namely that inconsistent reaction 
kinetics can contrive with skewed interaction Hamiltoni¬ 
ans to produce a deluding agreement with experiments. 

To our understanding, singlet and triplet projections 
are an integral part of the dynamics, the omission of 
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Kominis'trajectory Haberkorn's trajectory 



FIG. 3: (a)-(c) Quantum trajectories in Kominis’ approach, 
the initial state being |S) (8) I'fl')- (d)-(f) Quantum trajectories 
in Haberkorn’s approach, the initial state being |S)(8)|-IJ-). (a,d) 
singlet probability (Qs). Two projections to the singlet state 
are observed in this particular realization of Kominis’ trajec¬ 
tory. Such projections are missing from Haberkorn’s trajec¬ 
tories. (b,e) nuclear spin of the radical pair, After 

the instant of recombination, the trajectory contributes zero 
to all physical observables. (c,f) nuclear spin diz deposited to 
the reaction products at the instant of recombination. Since 
the considered radical-pair carries one spin-1/2 nucleus, the 
magnitude of diz correctly is 1/2. 


which is responsible for the inconsistent behavior of the 
conventional theory. The reason is that the Hamiltonian 
V coupling the spin degrees of freedom of the radical- 
pair to the vibrational reservoir (Fig. 1) is responsible 
for both effects, random projections and recombination. 
The former within 2''‘^-order perturbation theory and the 
latter through l®‘-order perturbation theory. As shown 
in [i^, both effects are of order V^. 

Clearly, there are still unresolved problems that must 
be addressed. While the quantum trajectory picture of 
our approach captures what we think are the underlying 
physics, our master equation does not do a perfect job 
in matching Monte Carlo. Hence more work is required 
to address this issue. Lastly, we recommend that the 
theoretical interpretation of a large number of CIDNP 
measurements performed over the last several decades 
should be meticulously revisited in light of a deeper un¬ 
derstanding of radical-pair reaction kinetics. 


Appendix A: Generation of quantum trajectories 

For each of the two initial states, we average 2 x 10® 
quantum trajectories. A single quantum trajectory is 
generated as follows. We split the time from t = 0 to 


^max = 20/fcT, at which point the reaction is practically 
over, into n = 5 x 10® steps of duration dt = tmax/n. 
We start with the initial state |'(/’o) at time f = 0, and in 
each time step we draw a random number x uniformly 
distributed between 0 and 1. In our approach we calcu¬ 
late the probabilities dps, dpT, drg and drx, and split 
the real interval [0,1] into five intervals. If the random 
number x falls within the 

• F* interval of length dps, we realize K1 and move 
on to the next time step 

• 2'^'® interval of length dpT, we realize K2 and move 
on to the next time step 

• interval of length drs, we realize K3 and termi¬ 
nate the particular trajectory 

• 4*^ interval of length dr^, we realize K4 and termi¬ 
nate the particular trajectory 

• 5*^ interval of length 1 — dps — dpT — drs — dr^, we 
realize K5 and move on to the next time step. 

In Haberkorn’s approach we calculate the probabilities 
drs and drx and split [0,1] into three intervals. If the 
random number x falls within the 

• 1®* interval of length drs, we realize HI and termi¬ 
nate the particular trajectory 

• 2“'® interval of length drx, we realize H2 and ter¬ 
minate the particular trajectory 

• 3’'*^ interval of length 1 — drs — dr^, we realize H3 
and move on to the next time step 

Examples of single quantum trajectories are shown in 
Fig 13] For vizualizing the dynamics, we also plot the evo¬ 
lution of the singlet probability, (Qs), depicting S-T os¬ 
cillations. Along those oscillations there are random sin¬ 
glet/triplet projections in Kominis’ trajectories, whereas 
they are completely missing from Haberkorn’s trajecto¬ 
ries. In both theories, the trajectory is terminated by the 
recombination event. 


Appendix B: Accuracy of calculations 

The accuracy of the master equation results (first line 
of Table HI) is limited just by the number of time steps 
n. By increasing n beyond 5 x 10®, the results converge 
to numbers within 1% of the ones stated here. 

In order for the code to run in a practical amount of 
time, we limited the Monte Carlo simulation (second line 
of Table HI) to 2x 10® trajectories. The whole simulation 
(the four Monte Carlo traces of Fig. 2) takes about 9 h in 
a 4-core machine running at 2.5 GHz. The running time 
of propagating the density matrix discussed previously 
is negligible compared to one Monte Carlo trace, since 
the former involves just one time propagation of an 8 x 8 
matrix, compared to 2 x 10® time propagations of an 8 x 1 
vector required for the latter. 
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The accuracy of the integral f diz is easily estimated 
by taking e.g. 100 consecutive points along a relatively 
flat part of the traces in FiglJl plotting the y-axis valnes 
in a histogram, and fitting with a gaussian. The relative 
error is at the level of 10%. Thus the positive sign of 
Haberkorn’s Monte Carlo result is correct to within lOcr. 
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